Mott Gap Excitations and Resonant Inelastic X-Ray Scattering in Doped Cuprates 
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Predictions are made for the momentum- and carrier-dependent degradation of the Mott gap 
upon doping in high-T c cuprates as would be observed in Cu K-edge resonant inelastic x-ray scat- 
tering (RIXS). The two-dimensional Hubbard model with second- and third- nearest-neighbor hop- 
ping terms has been studied by numerical exact diagonalization. Special emphasis is placed on 
the particle-hole asymmetry of the Mott gap excitations. We argue that the Mott gap excita- 
tions observed by RIXS are significantly influenced by the interaction between charge carriers and 
antiferromagnetic correlations. 
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Asymmetry in the electronic properties between 
electron- and hole-doped high-Tc cuprates is one of the 
key issues relevant to an understanding of the mechanism 
of high- T c superconductivity. In the cuprates, the antifer- 
romagnetic (AF) correlations are stronger in the electron- 
doped than hole-doped systems Q : In the electron-doped 
Nd2- rr Ce a; Cu04 (NCCO), the antiferromagnetism con- 
tinues up to x — 0.15, whereas it disappears for a much 
smaller x (~ 0.02) in hole-doped L^-^Sr^CuC^ (LSCO). 
The kind of carriers is reflected not only in the mag- 
netic properties but also in the single-particle excitations. 
Angle-resolved photoemission spectroscopy (ARPES) ex- 
periments revealed that while the minimum energy ex- 
citation appears at around (7r/2,7r/2) for hole doping, 
it is seen near (7r,0) for electron doping 2J. On the 
other hand, the charge-transfer excitations from the oc- 
cupied Zhang-Rice singlet band (ZRB) Q composed of 
Cu 3d x 2_ y 2 and O 2p a orbitals to the unoccupied up- 
per Hubbard band (UHB), i.e. Mott gap excitation, seen 
in the optical conductivity, shows similar behavior for 
electron and hole dopings In the optical process, 

the excitations occur without the momentum transfer to 
the electron system. Resonant inelastic X-ray scattering 
(RIXS) is now gaining importance as a powerful tech- 
nique for the investigation of the momentum-dependent 
excitations with the energies of the order of the Mott 
gap energy d SHU Efl E3. In particular, the Cu 
if -edge RIXS measurements on the parent materials of 
high-Tc cuprates [(J, revealed the characteristic mo- 
mentum dependence of the Mott gap excitation [ll| . It 
is, thus, interesting to clarify whether a doping asymme- 
try also exists in this high-energy region and how the na- 
ture of the Mott gap excitations in doped cuprates looks 
in RIXS. 

In this Letter, we demonstrate theoretically the differ- 
ence of the Cu K-edge RIXS between hole- and electron- 
doped cuprates. The single-band Hubbard model with 
long-range hoppings is adopted to describe the ZRB and 
UHB by mapping ZRB onto the lower Hubbard band 
(LHB) in the model. Then, the Cu Is and 4p orbitals are 
incorporated with the model to include the ls-core hole 
and excited 4p electron into the intermediate state of the 
RIXS process. The RIXS spectra are calculated by using 



the numerically exact diagonalization technique. We find 
that the excitation spectrum from the LHB to the UHB 
becomes broad and less momentum dependent upon hole 
doping. This is in contrast to the electron-doped case, 
where the momentum dependence of the spectrum of un- 
doped system remains, except that along the (1,0) di- 
rection. The difference in the spectra between hole- and 
electron-doped systems follows the carrier-dependence of 
short-range AF spin correlation. We also find that the 
momentum dependence along the (1, 0) direction in elec- 
tron doping comes from the fact that the doped electrons 
occupy the states at around k = (n, 0) in the upper Hub- 
bard band. 

By mapping the ZRB onto the LHB, which is equiva- 
lent to the elimination of O 2p orbitals [ll| , the Hubbard 
Hamiltonian with second- and third-nearest-neighbor 
hoppings is written as, 

H 3d = -* E 4<a*-*' E 4<a* 

(ij)l8t,CT <ij>2nd,0' 

E d\ a d la +h.e. + uY,<AlM) 

(i,j>3rd,0- 1 

where d\ a is the creation operator of 3d electron with 
spin a at site i, nf a — d\ a di ta , the summations (i,j)i s t, 
(ij)2nd, and (i,j)3 r <j run over first, second, and third 
nearest-neighbor pairs, respectively, and the rest of the 
notation is standard. The on-site Coulomb energy U cor- 
responds to the charge transfer energy of cuprates. 

The schematic process of the RIXS for Cu K-edge is 
shown in Fig. n( a )- An absorption of an incident pho- 
ton with energy u>i, momentum Ki, and polarization 6i 
brings about the dipole transition of an electron from Cu 
Is to Ap orbital [process (i)]. In the intermediate states, 
3d electrons interact with a ls-core hole and a photoex- 
cited 4p electron via the Coulomb interactions so that the 
excitations in the 3d electron system are evolved [process 
(ii)]. The Fermi energies are located in the UHB and 
ZRB for electron- and hole-doped cases, respectively. As 
shown in Figs. IHb) and^Jc), the excitations occur within 
the band as well as between the bands. The latter corre- 
sponds to the Mott gap excitation. The Ap electron in the 




FIG. 1: Schematic picture of Cu if-edge RIXS process (a) 
and densities of states for electron- (b) and hole-doped (c) 
cases. The dipole transitions of an electron between Is and 
Ap orbitals occur with (i) absorption and (iii) emission of a 
photon. The 3d electrons are excited (ii) in the intermediate 
sate. The Fermi energies [denoted by dotted lines in (b) and 
(c)] are in UHB and ZRB (LHB) for electron- and hole-doped 
cases, respectively. 



intermediate state goes back to the Is orbital again and a 
photon with energy lu{, momentum Kf, and polarization 
6f is emitted [process (iii)]. The differences of energies 
and momenta between incident and emitted photons are 
transferred to the 3d electron system. 

In the intermediate states, there are a ls-core hole and 
a Ap electron, with which 3c? electrons interact. Since the 
ls-core hole is localized in a small radius of the Cu Is or- 
bital, the attractive interaction between the ls-core hole 
and 3d electrons is strong. The interaction is written as, 
-His-M = -(7 r yi; „ „, nf„nf„,, where nf „ is the number 

J- 3 ^ L ^1,(T,CT 1,CT 1,(7 ' 1,(7 

operator of ls-core hole with spin a at site i, and U c is 
taken to be positive. On the contrary, the Coulomb in- 
teractions related to the Ap electron are neglected since 
the Ap electron is delocalized ^lj. Furthermore, we as- 
sume that the photoexcited Ap electron enters into the 
bottom of the Ap band with momentum ko . Under these 
assumptions, the RIXS spectrum is expressed as, 



J(AK,Ao,) = ]T 



S k -K f 



xS(Auj - E a + Eq), 



(2) 



where H = H 3d + H ls _ 3d + H lsAp , H lsAp being composed 
of the energy separation £i s _4 P between the Is level and 
the bottom of the Ap band, AK = K; Kf , Aoj = u>\ — uif , 



s k a (Pkcr) 1S * ne creation operator of the ls-core hole 
(Ap electron) with momentum k and spin er, |0) is the 
ground state with energy Eo, \a) is the final state of the 
RIXS process with energy E a , and T is the inverse of 
the relaxation time in the intermediate state. The RIXS 
spectrum of Eq. is calculated on a 4 x 4-site cluster 
with periodic boundary conditions by using a modified 
version of the conjugate-gradient method together with 
the Lanczos technique. 
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FIG. 2: Resonant inelastic x-ray scattering spectra (RIXS) 
for Cu K-edge of doped Hubbard model with long-range hop- 
pings, (a) The electron doped case ((n) = 1.125). (b) The 
hole doped case ((n) = 0.875). The spectra of the elastic 
scattering process at AK = (0, 0) are not shown. The param- 
eters are U/t = 10, U c /t = 15, Y/t = 1, t'/t = -0.34, and 
t" /t = 0.23. The 8- functions (the vertical thin solid lines) are 
convoluted with Lorentzian broadening of 0.2t. Insets are the 
Cu Is absorption spectra with TxAs/i = Y/t = 1.0, and the 
incident photon energies uii's for RIXS are set to the values 
denoted by the arrows. 



The values of the parameters are as follows: t'/t = 
-0.34, t"/t = 0.23, U/t = 10, U c /t = 15, and Y/t = 1 
with t = 0.35 eV that are estimated from the analyses 
of ARPES data ^1 an d are the same as those used in 
Ref. [HI. 

Insets of Figs. GJa) and 12b) are the Cu Is x-ray ab- 
sorption spectrum (XAS) for electron- and hole-doped 
cases, respectively 13]. The spectrum is given by 



1. 



D(u) = -Im(0|s ko _K.<rPko.o-7r- 



1 



Eq — Ul — ITxAS 



X Pko,^ko-Kj°>> 



(3) 



where H is the same as that in Eq. (0. In the inset 
of Fig. |21 (a), there appear three peaks, i.e., peaks at 
around u> — £i s -4 P = — 26i, — 20t, and — 13t. The peaks 
at around — 20i and —13t are also seen in the undoped 
case The peak at around — 20i corresponds to a 

final state where the core hole is screened and thus the 
core-hole site is doubly occupied by 3c? electrons {U — 
2U C = — 20t). This final state promotes the excitations 
from LHB to UHB in the RIXS. The peak at around 
— 13£ represents unscreened core-hole state and mainly 
contains the configuration that a core hole is created at 
a singly occupied site (— U c = — lbt). The peak at around 
— 26i, which appears upon electron doping, corresponds 
to a final state where the core-hole is created at a doubly 
occupied site induced by electron doping (— 2U C = — 30i). 
In the hole-doped case in Fig. Hfb), a peak appears at 
around At, in addition to the peaks at — 20i and — 13i. 
The peak at around At mainly contains the configuration 
that the core-hole is created at an empty site induced by 
hole doping. Because we are interested in the excitations 
from LHB to UHB in the RIXS spectra, we set u>i's to 
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FIG. 3: Momentum dependence of the center of gravity of 
RIXS spectra associated with the excitation from LHB to 
UHB. The spectral weight is adopted in energy regions above 
At, 8t, and 6t for undoped ((n) = 1, denoted by circles), 
electron- (1.125, upward triangles), and hole-doped (0.875, 
downward triangles) cases, respectively. The right axis de- 
notes the relative energy from the energy at AK = (0, 0) 
with t = 0.35 eV. 



the energy of the peak at around — 20i as denoted by the 
arrows in the insets of Figs.[2a) and[3fb). 

Figure [2f a) shows the RIXS spectra in the electron- 
doped case ((n) — 1.125, where (n) is the electron con- 
centration per site). The spectra below Aui ~ St are as- 
sociated with the excitations within the UHB, and those 
above 8t are due to the excitations from LHB to UHB, 
i.e., Mott gap excitations. We find that, in the energy 
region above 8t, the spectra strongly depend on momen- 
tum, showing a feature that the spectral weight shifts to 
higher energy region with increasing |AK|. The feature 
is similar to that of undoped case [llj. We note that, 
when the incident energy oji is set to ~ — 26t which is the 
energy of the lowest-energy peak in XAS, the RIXS spec- 
tra appear only in the lower energy region below Aw ~ 8t 
(not shown) . This is because the excitation from a singly 
occupied state to a doubly occupied state is reduced in 
the intermediate state of the RIXS process. Therefore, 
by tuning u)i to a peak energy of U — 2U C ~ —20t, we can 
enhance strongly the RIXS spectra associated with the 
Mott gap excitations. Since there is no such an enhance- 
ment effect on the Mott gap excitation in the optical 
conductivity 0, 0| as well as electron-energy loss spec- 
troscopy, we can say that by utilizing the incident energy 
dependence, the RIXS is the unique technique to inves- 
tigate the Mott gap excitation even in doped materials. 

Figure [2Jb) shows the RIXS spectra in the hole-doped 
case ((n) = 0.875). The spectra above the energy ~ 6t 
are associated with the Mott gap excitations. Compared 
with the electron-doped case, the spectra show broad fea- 
tures: For example, at AK = (tt,tt), the spectrum ex- 
tends from Auj ~ 7t to ~ 14i, and the energy position of 
the maximum spectral weight (~ lOi) is lower than that 
of electron-doped case (~ 13i). The spectra at other 
AK's are also extended to wide energy region similar to 
that at (7r,7r), and the energy distributions of spectral 
weights seem to be rather independent of momentum. 



In order to compare the energy positions of the spec- 
tra, the momentum dependence of the center of gravity 
of the RIXS spectra associated with the Mott gap ex- 
citation is plotted in Fig. [21 where the spectral weight 
is adopted from the energy regions above At, 8t, and 6t 
for undoped ((n) = 1, denoted by circles) [ll|, electron- 
(1.125, upward triangles), and hole-doped (0.875, down- 
ward triangles) cases, respectively. We find that the en- 
ergy positions in the doped cases are shifted to the high- 
energy side compared with the undoped case. This is 
because, upon hole (electron) doping, the Fermi energy 
shifts to LHB (UHB) and the energy difference between 
states of the occupied LHB and unoccupied UHB be- 
comes large. In the undoped case, the center of grav- 
ity shifts to higher energy with increasing |AK|. It has 
been shown that this momentum dependence has good 
correspondence with the experimental data [9|- We find 
that the momentum dependence in the undoped case re- 
mains by electron doping, except along the (1,0) direc- 
tion where the energy difference between (7r/2,0) and 
(0, 0) is larger than that in the undoped case. The origin 
of the similarity of the momentum dependence will be 
discussed below. Here let us discuss the physics behind 
the behavior along the (1,0) direction. In the undoped 
case, the RIXS spectra have a characteristic feature along 
the (1,0) direction where the edge of the RIXS spectrum 
at AK = (7r/2,0) is rather lower in energy than that 
at (0,0) [ill: The edge of the RIXS at AK = (tt/2,0) 
comes from the excitation from the occupied k = (7r/2, 0) 
state to the unoccu pied k = (tt, 0) one which is the lowest 
in energy in UHB [ljj. This edge feature restrains the 
center of gravity at (7r/2,0) from shifting to the higher 
energy in the undoped case. Since the (tt, 0) state of UHB 
is occupied by electrons upon electron doping, the center 
of gravity at AK = (tt/2,0) in the RIXS shifts to the 
higher energy. 

In contrast, the momentum dependence becomes 
weaker by hole doping, as seen in Fig. 0. Let 
us make clear that the difference between electron and 
hole dopings comes from the lack of particle-hole sym- 
metry in cuprates but not from the RIXS process it- 
self in which the background configurations around the 
core-hole screened site are different between electron and 
hole dopings. On Fig. 0] is plotted the center of grav- 
ity of RIXS spectrum on the Hubbard model with only 
the nearest-neighbor hopping, which has the particle-hole 
symmetry. To obtain the RIXS spectra, the incident pho- 
ton energies u>i's are set to the peak energy of U — 2U C 
in XAS spectra for both electron- and hole-doped cases. 
We can obtain almost the same momentum dependence 
in both electron and hole dopings accompanied by the 
same energy shift from the undoped case, although the 
configurations in the intermediate states are different. 
Therefore, the spectra of the Mott gap excitations have 
almost the same features between electron and hole dop- 
ings when the model has the particle-hole symmetry. 

One possible origin of the fact that the momentum- 
dependent spectra remain upon electron doping, whereas 
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FIG. 4: The same as Fig.[Ubut for t' = t" = 0. The spectral 
weight is adopted in energy regions above (St, It, and It for 
undoped, electron-, and hole-doped cases, respectively. 

the broad feature appears in the hole-doped case, comes 
from the doping effect on the AF correlations. In the 
undoped case, it has been discussed that the AF correla- 
tion plays a crucial role in the RIXS excitation [lTl Il5| : 
The matrix elements of the excitation from the top of 
the LHB [k = (>r/2,7r/2)] to the bottom of the UHB 
[k = (ir, 0)] are almost zero due to the coherence factor 
arising from the AF order. This fact affects strongly the 
momentum dependence of RIXS in the undoped case. 
Upon doping, the AF correlation changes. However, it 
strongly depends on carriers in the presence of the long- 
range hoppings 0, 0, : With the parameters used 
in Fig. |2 the short-range AF correlation is kept in the 
electron-doped case, whereas the correlation is strongly 
suppressed in the hole-doped case. Thus, we can suppose 
that the RIXS spectra in the electron-doped system are 
similar to the undoped case, while the spectra in the hole- 
doped system are different. In fact, as discussed above, 
in the electron-doped case, the calculated RIXS spectra 
show the momentum dependence similar to that of the 
undoped case, but remarkably different in the hole-doped 



case. Therefore, it is natural to consider that the Mott 
gap excitation is significantly influenced by the magni- 
tude of AF correlations even in doped systems. 

Finally, we discuss the material dependence of RIXS. 
The parameters of the second- and third-nearest- 
neighbor hoppings, t' and t" , play a role in the material 
dependence of the electronic properties. The hopping 
parameters t' and t" in LSCO are smaller than those in 
Bi 2 Sr 2 CaCu 2 8 (BSCCO), Ca 2 Cu0 2 Cl 2 (CCOC), and 
NCCO 19]. Therefore the features of the RIXS spectra 
on BSCCO, CCOC and NCCO are expected to follow 
Fig. 20], while for doped LSCO, the feature may be 
rather similar to that in Fig. 0] J2l| . 

In summary, we have demonstrated theoretically the 
difference of the Cu K-edge RIXS between hole- and 
electron-doped cuprates by using the numerically exact 
diagonalization technique on small clusters of the single- 
band Hubbard model with long-range hoppings. We have 
found that, upon electron doping, the RIXS spectra along 
the (1,0) direction shift to the higher energy side than 
those in undoped case. In contrast to electron-doped 
case, the spectra for the excitations from the lower to 
upper Hubbard bands show the broad feature and less 
momentum dependence in the hole-doped case. The dif- 
ference in the spectra between hole- and electron-doped 
systems follows the carrier dependence of short-range AF 
spin correlation. The RIXS experiments in a variety of 
doped cuprates are desired |20t l2ll| . 
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